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Abstract 

The efficiency of macroscopic heat engines is restricted by the second law of 
thermodynamics. They can reach at most the efficiency of a Carnot engine. 
In contrast, heat currents in mesoscopic heat engines show fluctuations. Thus, 
there is a small probability that a mesoscopic heat engine exceeds Carnot’s 
maximum value during a short measurement time. We illustrate this effect using 
a quantum point contact as a heat engine. When a temperature difference is 
applied to a quantum point contact, the system may be utilized as a source 
of electrical power under steady state conditions. We first discuss the optimal 
working point of such a heat engine that maximizes the generated electrical 
power and subsequently calculate the statistics for deviations of the efficiency 
from its most likely value. We find that deviations surpassing the Carnot limit 
are possible, but unlikely. 


M. Biittiker was among the first scientists to realize that measurements of cur¬ 
rent fluctuations deliver most valuable information about the internal structure 
of mesoscopic conductors [l[. The measurement of shot noise 0 in a tunnel 
junction, for instance, may be used to determine the elementary charge of the 
charge carriers transferred through the circuit. Its measurement may as well 
serve to reveal the transmission probabilities of a multichannel mesoscopic point 
contact. 

The description of current fluctuations was later extended to full statistics of 
the charge transfer through a mesoscopic conductor Q. From an experimental 
point of view, current fluctuations are probably the easiest to measure. Nev¬ 
ertheless, statistics for a number of other mesoscopic physical quantities have 
also been investigated: among them, combined charge-phase statistics in the 
superconducting state 0 , waiting time statistics of a closed volume @ , voltage 
statistics on a current biased point contact 
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Recently, interest has shifted to energy transport through mesoscopic struc¬ 
tures. The study of energy transport is partially motivated by the possibil^ 
to use small circuits to convert local temperature differences into voltages [7|. 
Heat currents are subject to fluctuations as well |m. Such fluctuations have 
been theoretically studied in a number of situations TIE El [II El- Although 
the direct measurement of fluctuations in heat current is probably difficult, in¬ 
direct consequences of energy fluctuations have been observed experimentally 

13,111,113- 


In this paper, we apply the theory of heat current fluctuations to a question 
of rather conceptual than practical interest. We consider a mesoscopic heat 
engine that converts heat partially into electrical work. Since heat currents 
in mesoscopic devices fluctuate with time, any quantity derived from the heat 
currents will fluctuate as well. In particular, the efficiency of heat to work 
conversion will depend randomly on time and therefore may exceed the Carnot 
efficiency for a short time, not on average, but sometimes with a non-vanishing 
probability. It is the aim of this work to quantify this probability. We illustrate 
our discussion with a quantum point contact, a narrow constriction between two 
electrodes that shows quantized linear conductance. 


1. System and Formalism 


Our system of interest is a mesoscopic point contact coupling left (L) and right 
(i?) reservoirs that have in general different chemical potentials fxji and 
different temperatures T^, Tr. This is an experimentally relevant situation 
that has been the subject of recent works 


i3,[!l,[i3. 

will flow from the left to the right reservoir when Tl > Tr. The energy flow 


is accompanied by a charge flow against a difference in the chemical potentials, 
Ml < /ifl- This charge flow corresponds to the electrical work generated by the 
point contact. The inset of Fig. [T] shows the setup and the sign conventions. 
The heat extracted from the left reservoir El has a positive sign, the (smaller) 
heat evacuated into the right reservoir Eji has a negative sign. We take the 
generated work IT as a positive quantity: 


El Eh = IT > 0 . 


The efficiency of the heat engine is then defined as 

_ IT El + Eh 
El ■ 


( 1 ) 

( 2 ) 


We consider that both electrodes have well defined chemical potentials and 
temperatures. Hence, relaxation processes in the reservoirs are assumed to be 
fast compared to the time r that our measurement takes. At the end of the 
experiment we record the amount of heat El that is extracted from the hot 
reservoir and the amount of heat Eh that is dumped into the cold reservoir. Due 
to thermal and quantum fluctuations, the whole transfer process is probabilistic 
and described by a probability distribution Pt^El, Eh)- Often it is more useful 


2 












Figure 1: (Color online) Inset: A quantum point contact connected to two reservoirs at 
different temperatures Tl^ Tji and chemical potentials fiL, fiji may be used as a heat engine 
that converts heat partially into electrical work W. Main plot: Efficiency rj = WIE^ 
versus the applied temperature ratio T^jT^ at the optimal working point, compared to the 
ideal efficiency of a Carnot process and the Curzon-Ahlborn efficiency at maximum power. 


to use the corresponding cumulant generating function instead of 

P. Both quantities are linked by Fourier transformation 

Pr{EL,ER)r^ y . (3) 

The generating function describing the statistics of heat transfer in a general 
two-terminal conductor is given in Ref. Here, we adapt this generating 
function for our purpose: 

Sr{i^L,liR) = ^ y&ln{l + r/i(l - Jr) 

+ r(l - - l) } ■ (4) 

This expression applies to a spin-degenerate single-channel point contact with 
transparency F. It also contains the leads’ Fermi occupation factors Jl^r = 
1/[1 -I- For convenience, we hereafter choose units such that 

the Planck constant, the unit charge and the Boltzmann constant are equal to 
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one. Any cumulant of the distribution Pr{EL, Eji) can be obtained from the 
generating function by taking partial derivatives 


{{{E^riEnr)) 


gm+ng 




^L—iR—O 


(5) 


Thermoelectric effects necessarily require that the transmission probability 
r of the quantum point contact depends on energy e. A basic and convenient 
model for this dependence was proposed by M. Biittiker If the potential 
barrier creating the point contact is a saddle, the transmission probability of 
one single transmission channel reads 


r(e) 


1 

1 _|_ ’ 


( 6 ) 


where Eq is the potential at the saddle and ujz gives the energy width of the 
transition region. In this work we will use a simplified version. We assume a 
very long point contact, uj^ oo, and choose the energy scale such that Eq = 0. 
The transmission probability then jumps sharply from zero to one 



e < 0 
e > 0 


(7) 


when energy surpasses the threshold set by eg. 


2. Optimal Working Point 


For a given set of temperatures Tr and Tr we may optimize the chemical poten¬ 
tials fiR and fiR such that production of work is maximized on average. We call 
this optimized situation the working point of the point contact. The result is 
equivalent to that of the efficiency at maximum power which has been analyzed 
in detail in Ref. 2]| in the context of scattering theory of quantum transport. 
From Eq. o we have 


((W)) = {{Er)) + {{Er)), 


d{{w)) 


dfJ-LM 

with the heat extracted from the heat reservoir given by 


= 0 . 


pOO 

{{Er)) - r / dE{fR- /r) {e - hr) , 
^0 

and the generated work obeying the Joule expression 

poo 

((W)) ^t{hR- ml) / dE {Jr - Jr) , 
Jo 


( 8 ) 


(9) 


( 10 ) 


for the nonzero voltage difference jjR ~ Ml- Combining the derivatives with 
respect to both chemical potentials we find 


t!± 

Tr 


Mfl 

Tr 


( 11 ) 
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Since Tl > Tr and < fJ-R (charge current must flow against the potential 
from the hot to the cold reservoir), it follows immediately that both chemical 
potentials have to be negative (we set the Fermi energy Ep = 0). Substituting 
Eq. (ITT]) into ((TF)) and calculating the energy integral yields the intermediate 
result 

((IE)) ^ In (l + , (12) 

which still needs to be maximized with regard to We define x = and 

obtain a symmetric transcendental equation and its numerical solution 


(1 + x) ln(l + x) = —x ln(x), x ~ 0.318, fiL,R — ■ (13) 


Interestingly, both potentials lie below the subband energy Eq a-nd therefore 
the Fermi functions could be approximated as Boltzmann factors. This fact 
demonstrates that the operational principle of our device is based on a classical 
distribution of electrons. 

For a compact representation of the results, we introduce the numerical 
constants 

/»00 7 I 1 X 

?/=—ln(x) ln(l + x) ~ 0.316, z= —^^—0.612, - ~ 0.517. (14) 

z 

The mean work and extracted heat for the working point become 

{{w))^ = y- in - nf , ( 15 ) 

TT 

{{El))u, = in - n) [n + (i - f) tr] . (le) 

With these results we are ready to calculate the mean efficiency at maximum 
power extraction 


{{'n))w 


1 - 


Tr 

Tl 




0.517 


VC 

1.483 - 0.483r;c ' 


(17) 


in terms of the Carnot efficiency yc = 1 — Tp/Tp. Equation (ITTI) has been 


derived in Ref. 21] using nonlinear scattering theory. The mean efficiency 
reaches at most 50% of the Carnot efficiency. The efficiency of heat engines 
operated at maximum power is additionally limited by the Curzon-Ahlborn ef¬ 
ficiency [i^, VC A = 1 — a/Tr/Tl, which sets a stronger bound than the 

Carnot efficiency. An expansion of Eq. (ITT)) . iiv))w = 0.3497yc-|-0.11477p, shows 
that our result falls below the Curzon-Ahlborn efficiency and approaches t]ca 
only near equilibrium. 

For comparison, in Fig. [T] we show iiv))w (full line), rjc (dashed line) and 
rjcA (dot-dashed line) as a function of the temperature ratio. In the strong equi¬ 
librium case, iiv))w reaches roughly half the Carnot efficiency. We emphasize at 
this stage that we are considering an idealized model. The obtained efficiency 
serves as a reference point for the upcoming calculation. The experimentally 
measurable efficiency of a realistic constriction would be probably much lower, 
although quantum point contacts perform rather well as compared with other 
low-dimensional systems 


24 |. 
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3. Distribution of the Efficiency 

We have thus discussed the mean efficiency of average heat to average work con¬ 
version. However, the efficiency of a mesoscopic heat engine in general fluctuates 
because both work and heat fluctuate over time. Therefore, efficiencies larger 
than the Carnot limit are in principle possible [ 2 ^. To see this, we calculate 
the efficiency distribution by integrating the complete distribution Ej^) 

with respect to one variable 

Priv) - j cLElPAEl, -(1 - AEl) . (18) 

After representing the complete distribution by its generating function, Eq. Q , 
the energy integration can be carried out. This step fixes one auxiliary field 

= (1 — s-iid leaves us with a single integration 

PrA) - J . (19) 

We evaluate the last integral to exponential accuracy by deforming the integra¬ 
tion path in the complex plane and by using a saddle point approximation, i.e., 
replacing the integral by its value at the saddle point 

- mini5„ [S'^((l - AAr, iCfl)] ■ (20) 

The minimization is carried out numerically. Our results for the distribution 
PtA) are shown in Fig. [5] for different temperature ratios Tr/Tr. The chemi¬ 
cal potentials are taken such that the point contact is operated at its optimal 
working point as discussed in Sec. [H The inset of the figure displays the saddle 
point solutions ArA) corresponding to the different distributions. 

We observe in Fig.[2]that the probability peaks at the mean efficiency ((7?)),,,. 
At the larger corresponding Carnot efficiency, it shows a weak minimum |25l |. 
This minimum is most pronounced for strong nonequilibrium {Tr/Tr = 0.05). 
The Carnot limit is attained in the exponential tails of the distribution. It is 
therefore very unlikely that our device beats the second law. A read-off example 
illustrates that the expected effect is indeed weak. Let us consider a fictitious 
measurement lasting 10“® s and taking place at temperatures Tr = 1000 mK 
and Tr = 450 mK (purple dashed curve in Fig. [2]). For these values, the 
efficiency at maximum power given by Eq. (El) is {{v))w = 0.23 (the mean 
efficiency at this temperature ratio) whereas the Carnot efficiency is rjc = 0.55. 
Now, taking into account the units of the Boltzmann constant kg and Planck’s 
constant h, the probability Pr{l)dT] of an efficiency of 100% (larger than the 
Carnot value) is found to be roughly 10^^ times smaller than the maximal 
probability Pr{0.23)dri corresponding to {{'ri))w = 0.23. It would thus take a 
total measurement time of 1000 s to measure such a second-law violating event 
once. 
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Figure 2: (Color online) Main plot: Distribution Pt of the efficiency rj = WjEi^ on a 
logarithmic scale. The Carnot efficiency 1 — Tji/Ti^ can be exceeded with an exponentially 
small probability. Inset: Saddle point solution used to calculate the distributions shown 
in the main plot. 


4. Conclusions 


To sum up, we have considered a mesoscopic heat engine (a quantum point 
contact) coupled to hot and cold reservoirs and operating in a nonequilibrium 
steady state. We have determined the efficiency that maximizes the work extrac¬ 
tion in terms of the reservoir temperatures. Our result is based on a transmission 
function that depends on energy and is valid for noninteracting electrons. M. 
Biittiker soon recognized that in the nonlinear regime of transport electron- 
electron interactions should be taken into account within a gauge-invariant 
current-conserving scattering theory . When the potential landscape is self- 
consistently calculated, the transmision probabiliW becomes a function of the 
applied temperature and voltage differences [i^ 0, [23 ■ This effect has been 
neglected here, which is a valid approximation for quantum channels strongly 
coupled to a nearby gate [s^. Nevertheless, we leave the question open as to 
what extend our results will be modified in the presence of nonlinearities and 
interactions. 

Furthermore, we have discussed departures of the second law of thermody¬ 
namics due to fluctuations of the efficiency that exceed the Carnot limit. For a 
quantum point contact, we have determined the full distribution of the efficiency 
and have quantified the probability that a given measurement deviates from the 
mean efficiency value. Our results are thus relevant for a careful assessment of 
the heat-to-work conversion properties of quantum coherent heat engines. 
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